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A large spectroscopic survey is constructed of galaxies in the Coma cluster. 
The survey covers a wide area (1 deg 2 ) to deep magnitudes (R ~ 19.5), covering 
both the core (high density) and outskirts (intermediate to low density) of the 
cluster. The spectroscopic sample consists of a total of 1191 galaxies, of which, 
760 galaxies are confirmed members of the Coma cluster. A statistical technique 
is developed to correct the spectroscopic sample for incompletness. The corrected 
sample is then used to construct R-band luminosity function (LF) spanning a 
range of 7 magnitudes (—23 < M R — 5 log h e5 < —16) both at the core and 
outskirts of the cluster. The R-band LF for the entire Coma cluster, fitted to 
Schechter form, gives; M* R = -21.79^9 + 5 lo 9 h 65 and a = -1.18j£o£ 

Dependence of the LF on local environment in Coma is explored. The LFs are 
found to be the same, within the errors, between the inner and outer regions and 
close to those from recent measurements for field galaxies. This is remarkable 
given the variation in the spectral types of galaxies between field and cluster 
environments. The steep faint-end slope for the LFs, observed in previous studies 
using photometric surveys, is not found here. However, the LF in this study is 
only measured to Mr = —16, compared to much deeper limits (M# ~ — 12) 
achieved in photometric surveys. 



- 2 - 



The total B-band LF for the Coma cluster, fitted to a Schecter form is; = 
-19.95±J^ + 5 log h 65 and a = -0.96^;^- This also shows a dip at M B = -18 
mag., in agreement with previous studies. The implications of this feature are 
discussed. 

The LF is studied in B — R color intervals and shows a steep faint-end slope 
for red (B — R > 1.35) galaxies, both at the core and outskirts of the cluster. 
This population of low luminosity red galaxies has a higher surface density than 
the blue (B — R < 1.35) star-forming population and dominates the faint-end of 
the Coma cluster LF. 

It is found that relative number of high surface brightness galaxies is larger at 
the cluster core, implying destruction of low surface brightness galaxies in dense 
core environment. 

Subject headings: galaxies: clusters — galaxies: clusters: individual(Coma)— 
galaxies: distances and redshifts — galaxies:kinematics and dynamics 



1. Introduction 

Detailed knowledge of the luminosity function (hereafter LF), defined as the number 
density of galaxies with a given luminosity, is essential for any observational study of the 
formation and evolution of galaxies, and for constraining galaxy formation scenarios and 
models for large scale structure. Over the last three decades, extensive studies of LFs have 
been performed, both in clusters and general fields. Although the move from wide-area 
photographic plate surveys to high-performance, sensitive CCDs has greatly advanced the 
subject, there are still a number of questions unanswered, including: 

• Is a single parametric form for the LF a reasonable approximation over the entire 
luminosity range? 

• What is the dependence of the LF on environment and does this also depend on the 
selection characteristics of the sample (i.e. depth, completeness)? 

• How do the LFs for different Hubble types or color intervals compare? 



1 Based on observations made with the William Hcrschel Telescope operated on the island of La Palma 
by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias, and also the Anglo-Australian Telescope of the Anglo-Australian Observatory. 
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• Do the high and low surface brightness galaxies follow different LFs? 

• How do the LFs of giant and dwarf galaxies compare? 

To address these questions, one needs a wide-area survey, complete to deep magnitudes. 
Moreover, in comparing different LFs, one must minimise observational biases (i.e. different 
depths) and selection effects (i.e. incompleteness). These are some of the reasons for slow 
progress in this field. 

Study of the field LF requires a redshift survey, complete to some magnitude limit. Due 
to the small numbers of intrinsically faint galaxies in even the largest magnitude-limited 
redshift surveys, field LFs do not normally extend to faint magnitudes. While this can be 
avoided for cluster galaxies, as all the objects are roughly at the same distance, there are 
still problems at faint magnitudes due to contamination by background objects. Moreover, 
compared to isolated systems, cluster galaxies are subject to dynamical evolution, affecting 
the shape and characteristics of their LF. The former problem can be resolved by exploiting a 
deep spectroscopic sample in clusters (and so spectroscopically-confirmed cluster members), 
while the latter is best addressed by studying large areas around clusters, covering both low 
and high-density regions. 

In this paper, we aim to explore the above questions by constructing the LF in the 
Coma cluster, avoiding potential problems which have affected previous studies. This has 
become possible with the advent of large-format CCDs, allowing wide-area surveys to deep 
magnitudes. Combined with fiber-fed spectrographs, one could then construct statistical 
samples of spectroscopic data on galaxies to faint magnitudes. 

At a redshift of z — 0.023, Coma is the richest nearby cluster, providing a laboratory for 
studying evolution of different types of galaxies. Other nearby rich clusters, such as Virgo 
and Ursa Major, are relatively unevolved, not allowing a direct comparison of the properties 
of galaxies as a function of their local environment. A detailed study of Coma also provides 
a control sample to compare with more distant clusters. This improves previous works in 
many aspects, including: 

• wide- area CCD coverage of the cluster, extending to >ldeg. from its center (this also 
covers the NGC4839 group); 

• full spectroscopic information for galaxies in the sample, allowing a measure of the LF 
independent of uncertain background corrections; 

• inclusion of the faintest galaxies for which spectroscopic data can be obtained (R~19.5), 
allowing strong constraints on the shape of the LF at faint magnitudes; 
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• availability of accurate CCD photometry in both B and R bands, providing color 
information for galaxies in the sample; 

• well-defined selection function for the spectroscopic sample. 

This is the fourth paper in a series studying the nature of giant and dwarf galaxies 
in the Coma cluster, using photometric (Komiyama et al. 2001; paper I) and spectroscopic 
(Mobasher et al. 2001; paper II) observations. Analyses of the diagnostic line indices (Pog- 
gianti et al. 2001; paper III), the radial dependence of galaxy properties (Carter et al. 2002; 
paper V) and the dynamics of giant and dwarf populations (Edwards et al. 2002) are already 
performed. Future papers will present a spectroscopic comparison with intermediate red- 
shift clusters (Pogiantti et al. 2002; in prep.) and a study of the scatter and environmental 
dependence of the color-magnitude relation (Mobasher et al. 2002; in prep.). 

In the next section, observations and spectroscopic sample selection are discussed. The 
selection functions are derived in section 3. Section 4 presents the LFs and their dependence 
on different physical parameters. The results are then compared with other studies in section 
5, followed by a discussion of the results in section 6. The conclusions are listed in section 7. 
We assume a distance modulus of DM = 35.13 mag. for the Coma cluster, corresponding 
to H = 65 km/sec/Mpc. 

2. Observations and Sample Selection 

To address the aims of this study, photometric and spectroscopic observations were de- 
signed to survey a wide area of the cluster (consisting of both core and outskirt fields) to deep 
flux levels. This is then complemented by a sample of brighter galaxies with spectroscopic 
data, compiled from different sources. The result is a sufficiently large spectroscopic sample 
to allow a detailed study of the LF of galaxies. The procedure for selecting the spectroscopic 
sample and its properties are summarised as follows. 

A wide-area photometric survey was carried out by the MCCD (a mosaic camera con- 
sisting of 5 x 8 CCDs and a field of view of 0.5 deg 2 ; Sekiguchi et al. 1998) on the WHT. 
A total of five fields in the Coma cluster were surveyed, each 32.5 x 50.8arcmin 2 . The 
photometric surveys were performed in both B and R bands, and are complete to R = 21. 
Details of the photometric observations, data reductions, source extraction and the construc- 
tion of the photometric catalogs are given in Komiyama et al. (2001; paper I). Follow-up 
medium resolution (6-9A) spectroscopic observations, using WYFFOS on the WHT, were 
then performed on two fields of this survey: Comal (at the center) and Coma3 (south-west 
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of Comal, including the NGC4839 group). The spectroscopic sample has a magnitude limit 
of i?=19.75, with the spectra having large enough S/N ratios to allow accurate measurement 
of diagnostic line indices (Poggianti et al. 2001; Paper III). The coordinates of the center of 
fields and the total number of objects in the spectroscopic sample in each field is listed in 
Table 1. The spectroscopic sample selection, spectroscopic observations and data reduction 
are explained in detail in Mobasher et al. (2001; paper II), where the spectroscopic catalog is 
also presented. The criteria for selecting the spectroscopic sample in the central (Comal) and 
outskirt (Coma3) fields are the same, with no radially-dependent selection biases present. 
This sample will be refered to as the Deep Spectroscopic Survey (DSS). 

The DSS is complemented by another sample with a brighter spectroscopic magnitude 
limit (R ~ 18), based on the compilation from Colless & Dunn (1996) and another spec- 
troscopic survey of Coma galaxies, carried out using the 2-degree Field (2dF) spectroscopic 
facility at the AAT (Edwards et al. 2002). This sample is refered to as the Shallow Spectro- 
scopic Survey (SSS). The total number of galaxies in each field for the SSS sample is also 
listed in Table 1. Line indices are not available for the SSS sample. 

The spatial distribution of galaxies in the DSS and SSS are shown in Figure 1, with 
the location and size of the MCCD survey areas overlaid: the Comal field at the center of 
the cluster, the Coma2 field to the west of Comal, and the Coma3 field to the south-west 
of Comal. Comparison between the redshifts of galaxies common to both the DSS and SSS 
samples shows a mean difference of 28±9kms _1 between the estimated velocities in the two 
surveys. 

The final spectroscopic catalog is compiled by combining galaxies with available redshifts 
from the DSS and SSS samples, the B and i?-band CCD data (i.e. magnitudes and surface 
brightness values) from the MCCD photometric survey, and the line indices for galaxies in 
the Comal and Coma3 fields. For galaxies common to both spectroscopic surveys, the mean 
redshifts were calculated. The B and i?-band magnitudes in this study are measured over a 
circular aperture of radius 3 times the Kron radius (see Paper II). This has the advantage 

Table 1: Coordinates of the center of each field and the number of sources with spectroscopic 
data in the SSS and DSS samples 



Comal 
Coma2 
Coma3 



Field 



R.A. Dec. 
J2000 

12 59 23.7 28 01 12.5 
12 57 07.5 28 01 12.5 
12 57 07.5 27 11 13.0 



DSS 
Members 
189 



90 



SSS 

Total Members 
302 282 
97 

188 102 



Total 
369 
184 
148 
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Fig. I. — The distribution of galaxies with spectroscopic data from DSS (solid points); 
SSS (open circles) and galaxies in common to both surveys (dotted circles). The three 
fields with available spectroscopic data are shown as rectangles: Comal (central), Coma2 
(west of Comal) and Coma3 (south-west of Comal). The total number of galaxies and the 
coordinates of the center of each field are listed in Table 1. 
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of scaling the photometric aperture in proportion to the size of galaxies and is therefore less 
susceptible to under-estimating the magnitudes for larger galaxies, as is often the case if a 
constant aperture is used for all the objects at the same distance. The Kron magnitudes 
used here are close to total (see paper I for details). The B — R colors are also measured over 
circular apertures of 3 times the Kron radius and hence, correspond closely to total colors. 

Galaxies in the range 4000 < cz < 10000 km s _1 are considered to be members of the 
Coma cluster (Colless & Dunn 1996). The total number of galaxies identified as cluster 
members in each of the three fields here, are also listed in Table 1. In case of Coma2 
field, only spectroscopic data from the SSS are available. For this reason, the depth of the 
spectroscopic survey in this field is shallower than that in the other two fields, with a different 
completeness function, as discussed in the next section. 

The photometric properties of galaxies in the spectroscopic sample are explored using 
their B — R colors and i?-band effective surface brightness (fi c s) distributions in Figures 2 
and 3 respectively. The effective surface brightness is defined as the mean surface brightness 
within the effective radius (radius containing half the total light of the galaxy, with the latter 
assumed to be the luminosity inside an aperure 3 times the Kron radius of that galaxy- paper 
I) of galaxy. The effect of seeing on the observed surface brightness was considered and found 
to be negligible. 

There are four galaxies with B— R > 2 which are confirmed members of the Coma cluster 
(2 galaxies in Comal and one galaxy in each of Coma2 and Coma3 fields). Such red galaxies 
were previously considered to be background objects (Conselice et al 2002), with none so far 
identified as a Coma cluster member. Considering the surface brightness distributions, we 
find that cluster members have a relatively brighter effective surface brightness compared to 
non-members (when corrected for surface brightness dimming and the K-correction due to 
redshift), with the tail of the surface brightness distribution for cluster members extending 
to brighter values. The median effective surface brightness changes from fi e ff = 20.50 
mag/arcsec 2 in the cluster to 21.50 mag/arcsec 2 for the field galaxies. While this was known 
from previous studies, there is also evidence for a monotonic decrease in surface brightness 
from core (Comal) to outskirts (Coma3). 



3. Cluster membership 

The spectroscopic surveys are incomplete, so we do not know for every individual galaxy 
whether or not it is a member of the cluster. Therefore, in order to compute the cluster LF, 
we assume that the spectroscopic sample is representative, in the sense that the fraction of 



- 8- 



140 
120 
100 
80 
60 
40 
20 
70 
60 
50 
40 
30 
20 
10 




Coma 1 




.Coma 2 



.Coma 3 



J i T L i U 



r l 4^4- 



0.4 0.8 1.2 1.6 
B-R 



2.4 0.4 0.8 1.2 1.6 2 2.4 

B-R 



2.8 



Fig. 2. — B — R color distributions in the three fields for the entire DSS+SSS spectroscopic 
sample (dashed line), for spectroscopically-confirmed cluster members (solid line), and for 
non- members (dotted line). 



- 9- 




Fig. 3. — i?-band surface brightness distributions in the three fields for the entire DSS+SSS 
spectroscopic sample (dashed line), for spectroscopically-confirmed cluster members (solid 
line), and for non- members (dotted line). 
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galaxies that are cluster members (which depends on both magnitude and position) is the 
same in the (incomplete) spectroscopic sample as in the (complete) photometric sample. If 
this assumption applies to both the spectroscopic surveys we are using, then they can be 
straightforwardly combined to give the cluster membership fraction over the full magnitude 
range of the two surveys. 

If N c (m, r|spec) and N t (m, r|spec) are, respectively, the number of cluster members 
(based on their redshifts) and the total number of galaxies within the combined spectroscopic 
sample in given ranges of apparent magnitude and cluster radius, then the probability of 
a randomly-selected galaxy in the combined spectroscopic sample being a member of the 
Coma cluster, P(m, r|spec), is 



I v iV c (m,r spec) 

P(m,r spec) = — — : -. 

7V t (m,r|spec) 

We simplify this by taking as our radial bins the ranges covered by the three fields of the 
survey, so that we have 



D f i \ iV C) i(m|spec) 

Pi m spec = -, 

iv t) j(m|spec) 

where i — 1, 2, 3 corresponds to the Comal, Coma2 and Coma3 fields. 

The cluster membership fractions as functions of apparent magnitude are presented in 
Figure 4 and listed in Table 2 for each of the three fields. The uncertainties in these fractions, 
assuming Poisson statistics, are estimated as 



dPi(m | spec) 1 1 \i/2 

Pi (m | spec) N c>i (m | spec) 7V M (ra|spec) 

where iV C) j(m|spec) is a binomial variable with its corresponding variance, as discussed in 
the next section. The errors corresponding to each magnitude interval in each field are also 
listed in Table 2. The membership fraction is determined down to R ~ 19.5 for the Comal 
and Coma3 fields, but only down to R— 17.75 for Coma2 (since only the SSS is available in 
this field). As expected, the Comal field, being centrally located, has a higher membership 
fraction at all magnitudes than the Coma2 and Coma3 fields. Coma3 has a comparable 
membership fraction to Coma2 at all magnitudes, even though it is further from the centre 
of the cluster; this reflects the presence of the NGC4839 group in this field. 



Table 2: Membership fractions in R-band magnitude intervals for the three fields. Errors are 
given in brackets. Photometric errors (0.03 mag) are significantly smaller than the bin sizes 
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Fig. 4. — Cluster membership fractions as a function of apparent magnitude for Comal 
(solid line), Coma2 (dotted line) and Coma3 (dashed line). 
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4. The cluster luminosity function 

Let N P} i(m) be the number of galaxies as a function of magnitude m and field i, in the 
full photometric sample. The cluster LF can then be estimated as 

N Pti (m) 

<M m ) = — -J — fi( m ) , I 1 ) 

where is the area of each field, N Pii (m)/Ai the surface density in the photometric catalogue, 
and fi(m) = Pi (m| spec) is the cluster membership fraction. 

The errors in the LF can be estimated by noting that <pi{m) is the product of two 
random variables: N Pii (m), which is a Poisson random variable, and N c ^(m\spec) , which is a 
binomial random variable (since it is the number of cluster members found in a fixed number 
of trials, N ti i(m\spec) , where the probability of a single trial yielding a cluster member is 
the true cluster membership fraction, /j(m)). The relative error in the LF is given by the 
quadrature sum of the relative errors in N Pii (m) and N c>i (m\spec) , so that 

var^m)] _ vai[N Pti (m)] + var [A^(m| spec)] 
<pi{m) 2 N Pji (m) 2 N C:i (m\spec) 2 

The variance in the Poisson variable is its expectation value, vax[N Pti (m)] = N Pti {m), while 
the variance in the binomial variable is var [N c>i {m\ spec)} = N tii (m\spec)) /j(m)[l — fi(m)]. 
The relative error in the LF is thus 

5 Mm) = var^m)] 1 / 2 = / 1 + 1 1 

4>i(m) (t>i{m) \N Pii (m) N Cyi (m\spec) N tii (m\spec) 

This expression reduces to N c ^{m\spec)~ l l 2 , if N tji (m\spec) = N Pti (m) (i.e. if the spectro- 
scopic sample is complete), and to JVp^ra) -1 / 2 , if N Cji (m\spec) = N tji (m\spec) (i.e. if the 
membership fraction is unity). 

The LFs for the Comal, Coma2 and Coma3 fields are listed in Table 3 and presented in 
Figure 5. Magnitudes are corrected for both Galactic absorption (Ar = 0.02) and redshift 
(K z = 0.02) dimming. The distance of the Coma cluster is then estimated as m — M = 
DM + A R + K z where DM = 35.13 mag. is its distance modulus (corresponding to H = 65 
km/sec/Mpc). At this distance, 1° = 1.77 Mpc (assuming q = 0.5), corresponding to an 
area of 1.619 Mpc 2 for each of the three fields surveyed in this study. The LFs span the 
range —22 < M R — 5 log h 65 < —16, except for the Coma2 field, which is only based on the 
shallower SSS sample and hence, only extends to M R = —18. 

It is clear from Figure 5 that the bright-end of the LFs are not well constrained from 
these three fields (filled circles). This leads to instability in parametric fits to LFs, as the 
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Table 3: The observed LFs and their associated errors for the three fields. LFs are in 
galaxies/ mag/ Mpc 2 . Galaxies from both the extended SSS (lines 1-9) and DSS+SSS (lines 
10-24) are used. A distance modulus of DM = 35.13 (H = 65 Km/sec/Mpc) is assumed 
for Coma. 
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characteristic magnitudes and faint-end slopes are known to be correlated (see the next 
section for details). To overcome this, we construct the bright-end of the LFs by combining 
the present survey (DSS+SSS) with an extension of the SSS sample, covering a larger area 
of the Coma cluster. This extended survey covers an area of 1° radius (corresponding to 
the field of view of 2dF), centered on the Coma core, and to the same depth as the SSS 
sample, substantially increasing the number of bright galaxies in the spectroscopic Coma 
survey. Here, we assume that the shape of the bright-end of the LFs is independent from 
their local environment, with the implication of this assumption discussed in section 4.1. 
However, photometry for galaxies in the extended SSS sample is only available in b and r 
bands, measured from photographic plates (Godwin et al 1983). To convert these to the 
photometric system used in this study, we use galaxies with available photometry in both 
systems and derive the empirical relation 

R = r + k 1 (b-r) + k 2 

where k\ = 0.101 ± 0.068 and k 2 = 0.143 ± 0.037 with rms = 0.218 mag. This relation is 
subsequently used to convert magnitudes in the extended SSS sample to those of the MCCD 
photometric system of the present study. 

The membership fraction is estimated for this sample, as discussed in the last section, 
and used to construct the LF over a 1-degree radius region about the cluster center. The 
LF, binned in 0.25 mag. intervals, is then normalised to those of Comal, Coma2 and Coma3 
fields, derived from the DSS and SSS samples, using galaxies with M R < —21 mag., and 
shown by open circles in Figure 5. This extends the bright-end of the LFs to M R = —23. 
By adopting the extended SSS sample here, we impose strong constraints on the bright-end 
of the LF for each field while, at the same time, retaining the main advantage of the present 
sample in extending the Coma LF to the faintest possible magnitudes. The estimated bright- 
end of the LF, derived from the extended sample and normalised to the LFs in each of the 
three fields, are also listed in Table 3. In the next section, parametric fits are carried out to 
the observed LFs over the entire magnitude range shown in Figure 5. 

4.1. Parametric fits to luminosity functions 

We consider a Schechter parametric form for the observed LFs, expressed as 

$ S (M) = (/)* X a+1 e- x 

where X = 1 o-°- 4 ( M - M *) ; with M* (characteristic magnitude), a (faint-end slope) and 0* 
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Fig. 5. — Observed LFs for Comal, Coma2 and Coma3 fields. The DSS+SSS sample for a 
given field is shown by filled circles and the extended SSS sample (normalized to DSS+SSS 
survey) by open circles. Schechter LF fits, as listed in Table 3, are also shown. 1,2 and 3a 
error contours for M£ and a are presented on the separate panel. 
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(normalisation) the parameters to be determined (Schechter 1976). 

Galaxies from DSS and SSS, combined with those from the extended SSS sample, are 
used in the fit (both the solid and open circles in Figure 5). As is clear from Table 3 and 
Figure 5, galaxies from the DSS+SSS and extended SSS samples overlap in the magnitude 
range —22.50 < M R < —21.00. For galaxies over this interval, the average of 4>(M) values 
are calculated in each magnitude bin and used to fit the LF. This allows the magnitude 
bins to be independent in the fitting process. Using only one of the DSS+SSS or extended 
SSS samples over their common magnitude interval, will not affect the results. The fits are 
carried out over the luminosity range —23 < Mr — 5 log < — 16 for Comal and Coma3 
and —23 < Mr — 5 log < —17.5 for Coma2 fields. The fits for individual fields are shown 
in Figure 5 and presented in Table 4. The reduced x 2 values, x 2 l v i number of degrees of 
freedom, v, and the goodness of fit probability, P(x 2 \ iy ), are also listed for each field. The 1,2 
and 3cr error contours corresponding to each of the LFs are shown in Figure 5, with the la 
error estimates for LF parameters presented in Table 4. The x 2 results in Table 4 indicate 
Pix 2 ^) > 20% probability that the data, in all the three fields, are drawn from a Schechter 
LF, confirming that a Schechter form provides acceptable fits to individual LFs both at the 
core and in the outskirts of the cluster. 

The procedure of using the bright-end of the LF from the entire Coma cluster to con- 
strain that of the individual fields (located at different positions in Coma), makes the implicit 
assumption that the LFs at bright magnitudes are independent of their local environment. 
To explore the extent to which this assumption could affect the overall shape of the LFs and 
their parametric form, we now fix M* values for each field to those measured and listed in 
Table 4, and carry out a two parameter fit to DSS+SSS (ie. the sample which extends to 
fainter magnitudes) to determine a and 0*. These results together with their corresponding 
goodness of fit estimates are also presented in Table 4. There is no difference (within the 
errors) in the faint-end slope and normalisation of LFs between the two cases, with both 
providing acceptable fits to the observed LFs. The probability that these data are drawn 
from a parent Schechter LF is P(x 2 |z/) > 10%, further confirming the previous result that 
a Schechter form provides acceptable fit to the LF at the core and outskirts of the Coma 
cluster. This justifies the above assumption that the bright-end of the LFs are similar among 
the three fields in the Coma. 

Using the techniques developed and the results so far, we now concentrate on addressing 
the following three questions: 

(i) Is the LF the same at the core and outskirts of the Coma ? how does it depend on 
luminosity? 
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Table 4: Results from the parametric fits to the observed LFs in Comal, Coma2, Coma3 and 
the entire Coma cluster, la error estimates are also listed. For each of the three fields, the LF 
parameters are estimated after constraining their bright-end using data from the extended 
SSS sample (1st line for each field). Then the Mjt, values are fixed to these estimates and a 
2-parameter fit is performed to measure a and <p* (2nd line for each field). ' Coma3 LF with 
the NGC4839 group removed. 
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(ii) Can the LF for the entire Coma cluster be parametrised by a single form over its 
entire luminosity range? 

(iii) Does the LF change, depending on the observed wavelength or color and surface 
brightness of galaxies? 

4.2. Environmental dependence 

While the characteristic magnitudes for the LFs in the Comal, Coma2 and Coma3 fields 
are close (Table 4), there is evidence for a slight increase in the faint-end slope from —1.16 
at the core (Comal) to —1.29 at the outskirts (Coma3). However, this is only a la effect. 
It is also clear from the x 2 tests that a Schechter function form provides acceptable fits to 
the observed LFs in all three fields and over the entire magnitude range covered. Moreover, 
this shows that the three LFs have the same shape over their common luminosity range. 

Given the similarity of LF shapes in different environments in the Coma cluster, it is 
appropriate to sum the individual LFs to estimate the total LF, $(M), for the Coma cluster 

as 

$(M) = £? =1 &(M) 

While this provides an accurate estimate of the faint-end (Mr > —18), it is only weakly 
constrained at the bright-end (Mr < —21), with large Poisson errors. Therefore, the total 
LF is derived here by combining galaxies in the DSS+SSS samples (Comal, Coma2, Coma3 
fields), as measured by $(M), with those in the extended SSS (all galaxies in the shallow 
spectroscopic survey covering an area of 1° diameter in the Coma). The membership fraction 
is measured as discussed in the last section. The two LFs are normalised over the magnitude 
range -22.5 < M R < -21, with only the extended SSS sample at M R < -17.75 used 
for fitting to the parametric form. The total LF is then derived over a magnitude range 
—23 < Mr — 5 log < — 16 and presented in Figure 6. The total LF over a 1-degree 
radius region about the cluster center is listed in Table 5, with its corresponding parametric 
fit presented in Table 4. 

The reduced x 2 value, corresponding to the LF from all fields (last line in Table 4), 
indicates a small probability (P(x 2 \v) ~ 1%) that the data are drawn from a Schechter LF. 
We now consider the possibilities which might contribute to this. The total LF in Figure 6 
is constructed assuming the LF for galaxies covering the entire 1° diameter in the Coma to 
have the same shape as those in the three (Comal, Coma2 and Coma3) fields. However, it is 
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likely that this is not the case for other regions of the Coma cluster (due to differences in local 
densities), affecting the bright-end of the LF. Indeed, constraining the sample to only the 
three fields studied here (DSS+SSS galaxies located in Comal, Coma2 and Coma3 fields), 
increases the probability to P(x 2 \v) ~ 10%, with no significant change in the estimated LF 
parameters listed in Table 4. Moreover, differences in normalisations ((f)* values) between 
the three fields are likely to result in a poor fit to a Schechter LF form once the fields are 
combined. Alternatively, it is possible that a 2-component LF (ie. Gaussian+Schechter) is 
needed to model the total LF for the Coma cluster (Biviano et al 1995; Yagi et al. 2002). 

The observed faint-end slope of — 1.18l 'g 2 f° un d here, while consistent with individual 
fields (Comal, Coma2, Coma3), is relatively shallower than other studies of Coma cluster LF 
(Trentham 1998; Seeker & Harris 1996; Beijersbergen et al. 2002). However, the spectroscopic 
LF here has a much brighter magnitude limit (Mr = —16). Moreover, the measurements 
based on photometric surveys are mostly unconstrained at faint magnitudes due to small 
number statistics, incompleteness, uncertain background correction and contamination by 
globular clusters (see section 5 for details). 

Study of the dependence of LFs on their local environments here could be hampered 
by the presence of the NGC 4839 group in the Coma3 field. This is examined by removing 
galaxies associated with this group from the Coma3 sample and then, measuring its LF. 
Defining NGC 4839 group as galaxies within a circle of radius 9 arcmin. centered on the cD 
galaxy NGC 4839 (Beijersbergen et al. 2001), we construct the Coma3 LF after removing 
these objects and then carry out a parametric fit to the Schechter LF form. The result is 
also listed in Table 4 and shows a slight increase (la) in the faint-end slope of the Coma 
LF from the core to outskirts. Spectroscopic data on another field in the outskirts of the 
Coma cluster at a similar distance from the center as Coma3 is needed to accurately measure 
changes in the faint-end slope with radius. 

4.3. Wavelength Dependence 

The B-band LF for the entire Coma cluster is constructed and presented in Figure 7, 
following the same procedure discussed in section 4. The extended SSS (covering 1 deg 2 
area of the Coma) and DSS samples are used with the completness functions estimated in 
the same way as for the R-band LF. The B-band magnitudes which were not available for a 
fraction of galaxies in the extended SSS survey are estimated using their R-band magnitudes 
(measured from photographic plate b and r— band data in section 4) and B — R colors. The 
total B-band LF is also tabulated in Table 5. 
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Table 5: Total LF (<E>(M)) for the entire Coma cluster. Galaxies from the SSS sample (Mr < 
-17.5; M B < -16.25) are combined with data from DSS (M R > -17.5; M B > -16.25). 
Errors are estimated as discussed in the text. The last four lines correspond to the DSS 
survey. 
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Fig. 6. — Total LF, covering an area of 1° radius of the Coma cluster. Combined sample from 
the SSS (open circles) and DSS (filled circles) are used. Schechter LF fit is also presented, 
and a error countours at 1,2 and 3<j levels are also shown. 
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The spectroscopic sample in this study was selected in R-band. This was motivated 
because the light at longer wavelengths is dominated by evolved stellar population and 
because the Coma cluster is rich in early-type evolved galaxies. Therefore, we expect a bias 
against the very blue and faint galaxies, leading to a flatter faint-end slope for the B-band LF 
here. Having this in mind, the B-band LF is fitted to a Schechter form, with the best fit and 
1,2 and 3a error contours presented in Figure 7. We estimate M B = —19. 951^10 + 5 ^655 
a = -0.96l[};{£; 0* = 16.3±g;i h 2 m Mpc' 2 and ^jv = 1.16. This corresponds to P(x» = 
25% probability that a Schechter LF form is a good representation of the data. 

The total B-band LF in Figure 7 shows a dip at M B = —18 mag. A similar gap 
with smaller amplitude has been observed in the total R-band LF (Figure 6) at almost 
the same magnitude, = —19.5, shifted by the mean color of galaxies in this sample 
(< B — R > = 1.5). Although this is only a la effect and likely caused by Poisson statistics, 
a detailed study of this feature is useful as it gives clues towards the overall shape of the LF. 
Such a behavior has also been confirmed from other studies at almost the same magnitude 
in the rest-frame B-band LFs of both nearby (Biviano et al 1995) and intermediate redshift 
(Dahlen et al 2003) clusters. Studying the R-band LFs for a photometric sample of 10 clusters 
at different redshifts and with different richness classes, Yagi et al (2002) show evidence for 
a dip in the cluster LFs at the same magnitude as here. After dividing their sample into 
type specific groups, consisting of elliptical (r 1 / 4 -like profile) and spiral (exponential-like 
profile), they conclude that the observed dip is almost entirely due to contribution from 
the early-type (r 1//4 -like) galaxies to the total LF. They also find that clusters with larger 
velocity dispersions have more distinct dips, perhaps indicating that the amplitude of the dip 
depends on the dynamical state of the cluster and consequently, on the dominant population 
of galaxies in that cluster. 

To explore the effect of this gap on the general shape of the LF, we now remove galaxies 
in the range —18.5 < M B < —17.75, which contribute to the gap in Figure 7 and then fit 
the LF to a Schechter form. We find M* B = -19.74 ±0.10 + 5 log h 65 ; a = -0.90 ± 0.02; 
(f)* = 20.7 ± 0.1 h 2 Mpc~ 2 and x 2 l v — 0.97. This increases the probability of a Schechter 
LF form being a good representation of the data to P(x 2 ] v) — 50%. Given this result, it is 
likely that a two-component shape for the LF, consisting of a Gaussian, representing giant 
galaxies (M B < —18), and a power-law, representing the dwarf population (M B > —18), 
may well be appropriate for the composite LF over its entire magnitude range. 



-26- 




Fig. 7. — Total B-band LF, covering an area of 1° radius of the Coma cluster. Combined 
sample from the SSS (open circles) and DSS (filled circles) are used. Schechter LF fit is also 
presented. and a error countours at 1,2 and 3a levels are also shown. 
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4.4. Color Dependence 

The survey performed here is extensive enough to allow a study of the LFs both in 
different environments in the Coma cluster and in color intervals. This requires knowledge 
of the membership fraction in B — R color intervals. To minimise statistical uncertainties, we 
use galaxies from all the fields (in both DSS and SSS samples) to derive a common selection 
function, assuming no spatially-dependent selection bias. 

Two color intervals are adopted, corresponding to blue (B — R < 1.35) and red (B — R > 
1.35) galaxies. This color is chosen so that it is consistent with the average color of an 
intermediate type spiral (ie. Sc), allowing for clear separation between star-forming (blue) 
and evolved (red) populations. Moreover, this provides sufficiently large number of blue and 
red galaxies for statistical analysis. The blue galaxies here are also found to have strong 
emission lines, further confirming that they are undergoing star formation activity. The 
membership fractions, derived for the two color bins following the procedure explained in 
section 3, are presented in Figure 8. It is clear that, for both blue and red galaxies, the survey 
is 80% complete to R ~ 16, droping to 40% completeness at R ~ 18.5. Using membership 
fractions, the LFs are derived in B — R < 1.35 (star- forming) and B — R > 1.35 (evolved) 
intervals and presented in Figure 9 for the three fields separately and all the fields together. 

The evolved galaxies in Comal have a surface density which is an order of magnitude 
higher than that of star-forming galaxies, with the difference decreasing towards the outskirts 
(Coma2 and Coma3 fields). The sudden increase in the slope of the Coma3 field LF for red 
galaxies at M R > —18 is due to the contribution from the NGC4839 group. The B — R 
color distribution for spectroscopically confirmed members of this group (Figure 10) shows 
a significant number of these galaxies having B — R > 1.35. The similarity between the 
slope of the LFs for star-forming galaxies in different environments contradicts the result 
from Beijersbergen et al. (2002) who found a relatively steeper U-band LF in the outerpart 
of the Coma cluster. The difference here is likely due to selection of the present sample in 
the red passband, statistical fluctuations in both studies and uncertainties in background 
subtraction in the Beijersbergen et al. sample. 

Recently, it was shown that for a given spectral type of galaxies, the faint-end slope of 
the cluster LF is steeper than that in the field (De Propris et al 2003). This is in apparent 
contradiction with the result here, in which similar LFs are found for both the blue and red 
populations in the core (Comal) and outskirts (Coma3) of the Coma. However, clusters 
used in De Propris et al (2003), cover a wide range in richness, with many being as rich as 
the Coma3 field here. This indicates a larger density contrast between the cluster sample of 
De Propris et al and the general field, compared to the Comal and Coma3 fields. Moreover, 
the results from the two studies are not directly comparable as there is only a loose relation 
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Fig. 8. — Membership fractions in B — R color intervals for galaxies in DSS+SSS samples; 
B - R > 1.35 (solid line) and B - R < 1.35 (dashed line). 
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Fig. 9. — LFs in B — R color intervals for galaxies at the core and outskirts of the Coma 
cluster. B — R > 1.35 (filled circles) and B — R < 1.35 (open circles). The errorbars are 
estimated as discussed in the text. 
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between colors and spectral types used in these two studies. 

4.5. Surface Brightness Dependence 

The surface brightness distributions in Figure 3 show that, on average, the cluster mem- 
bers have a relatively brighter surface brightness compared to field galaxies. This implies 
that higher surface brightness galaxies tend to reside in environments with relatively larger 
densities. Figure 11 presents luminosity distributions for spectroscopically confirmed mem- 
bers of the Coma cluster in three surface brightness intervals; // e // < 20, 20 < /i e // < 22 
and fi e ff > 22 mag./arcsec 2 . There is a clear trend in the sense that low surface bright- 
ness galaxies have fainter magnitudes, reflecting a monotonic magnitude-surface brightness 
relation. Such relation has also been found for field galaxies (Brown et al 2001; Cross et al 
2001). However, there is also a considerable number of low surface brightness (/i e // > 22) 
galaxies in all the Coma fields (> 50%)- (Figure 11). This indicates a low surface brightness 
population that dominates the faint-end of the LF, in agreement with results from Andreon 
& Cuillandre (2001) and Sprayberry et al (1997). 

5. Comparison with other Studies 

Recent studies of the Coma cluster LF have led to steep faint-end slopes in both optical 
(Thompson and Gregory 1993; Seeker & Harris 1996; Trentham 1998) and near-IR (Mobasher 
& Trentham 1998; De Propris et al. 1998; Rines et al. 2001) wavelengths. All these studies are 
based on photometric surveys with statistical background corrections from control fields and 
varying fields of view. Recently, this technique has gained popularity by the advent of large- 
format and highly sensitive panoramic detectors, allowing more statistically representative 
measure of background contamination. 

Using photographic plates, Thompson and Gregory (1993) determined the LF to B ~ 20 
(Mb ~ —15). They found a steepening of the slope at faint magnitudes (a = —1.4). This 
was followed by Bernstein et al. (1995) study, who used CCD surveys, extending the Coma 
LF to much fainter limits (Mr ~ —9.5) over a smaller region. These authors also found a 
steep faint-end slope (a = —1.8). However, their result was significantly affected by globular 
cluster contamination in the halo of NGC4874, making their LF unconstrained fainter than 
M R = — 14. A deeper survey over a much larger area (700 arcmin 2 ) by Seeker & Harris 
(1996) also found a LF that rises to M R = — 14 and flattens thereafter. The flat bright-end, 
followed by steep faint-end slope was further confirmed by Lobo et al. (1997). 
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Fig. 10. — B — R color distribution for spectroscopically confirmed members of the NGC4839 
group in the Coma3 field. 
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Fig. 11. — Absolute magnitude distributions in effective surface brightness intervals for 
galaxies at the core, outskirts and the entire Coma cluster; fj, e ff < 20 (solid line), 20 < 
fx e ff < 22 (dotted line) and /i e ff > 22 (dashed line) mag/arcsec 2 . 
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Recently, Beijersbergen et al. (2002) carried out a multi-wavelength (UBr) study of LF 
in the Coma cluster. They found M R = — 21.47+q;]; 2 . + 5 log h 65 and a = —1.16, with a 
marginally steeper faint-end slope at larger radii from the cluster center. The increase in 
the faint-end slope of U-band LF was interpreted as due to a population of star-forming 
dwarf galaxies. No such increase in the number density of this population is found, as is 
clear from the B-band LF in Figure 7. Also, Andreon and Cuillandre (2001) carried out 
a wide-area (650 arcmin 2 ) photometric survey of the Coma cluster. Combined with HST 
observations, they confirmed that unresolved globular clusters could be mistakenly classified 
as dwarf galaxies at faint magnitudes, causing the observed steepness of the LF faint-end 
slope. Correcting for this, they derived a = —1.4, with both the shape and slope of the LF 
depending on color. They also derived LFs in surface brightness intervals and found low 
surface brightness galaxies as the main contributor at fainter magnitudes, as expected from, 
for example, studies of the local group (Mateo 1998). This result is in close agreement with 
that found in section 4.5, using a spectroscopic sample, and that in Trentham (1998) from 
a photometric survey. 

Trentham (1998) avoided some of the problems affecting the above studies by using a 
deeper control field to improve background estimates and contamination (at the faint-end) 
by globular clusters. Surveying a large area (674 arcmin 2 ) to deep (R = 24) limits, he 
found a steep faint-end slope (a = —1.7) at R > —17.6 followed by a turn over in the LF at 
R = —15.6, in both the core (r < 200 kpc) and outskirt (r > 200 kpc) fields. The population 
of galaxies at R ~ 19 were found to be dwarf spheroidals with B — R — 1.3. The narrow 
spread in B — R colors here was taken as evidence for a homogeneous population of dwarfs 
in Coma, in agreement with recent results. However, the slope and general shape of the LF 
at faint magnitudes only depend on a few points with large errorbars, making this part of 
the LF essentially unconstrained. 

All the above studies are based on photometric surveys with statistical background 
corrections. The faint-end slopes of the LFs from photometric surveys discussed above, are 
relatively steeper than that found from the present study, based on a spectroscopic survey. 
This is mostly due to the depth of the photometric surveys, extending to Mr ~ —12, as 
compared to the shallower spectroscopic survey here (Mr < —16), with the up-turn in the 
faint-end slope appearing at M R > —16. Moreover, it is possible that a combination of 
uncertain background correction and contamination by globular clusters, specially at faint 
magnitudes, may have led to the steep faint-end slopes found in studies based on photometric 
surveys. However, in a recent study, using deep simulated catalogues of cluster galaxies 
constructed from a flat LF (a = —1), Valotto et al (2001) show strong tendency towards 
steep faint-end slope for LFs (a > —1.5) in clusters selected in two dimensions. This indicates 
the cause of the observed steep faint-end slope in clusters, seen in photometric surveys, as 
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due to projection effect. 

Using a sample of 205 spectroscopically confirmed galaxies with 6 2 6.5 < 18 at the 48 x 25 
arcmin 2 center of the Coma cluster, Biviano et al (1995) found that a single Schechter 
function cannot adequately fit the LF of the Coma cluster. He found that a combination of 
Schechter and Gaussian forms provides a significantly better fit to the LF. However, Biviano 
et al's data do not strongly constrain the faint-end of the LF, with his spectroscopic sample 
becoming progressively incomplete towards fainter magnitudes. 

It is useful to compare the Coma LF derived here, with those in other rich clusters 
(Virgo), less evolved clusters (Ursa Major) and the general field. A recent study of the 
Virgo LF finds a steep faint-end slope of a = —1.6, based on photometric B-band CCD 
survey (Trentham & Hodgkin 2002). Although steep, this is significantly shallower than the 
slope of a — —2.2, found in i?-band by Phillipps et al. (1998). Despite the fact that these 
are based on power-law fits over restricted magnitude range, both the studies show the LF 
becoming relatively flat after the rise at the bright-end and before steepening towards fainter 
magnitudes. However, due to dominance of the star over galaxy counts and steep galaxy 
counts in the Virgo, measurement of the LF for this cluster is uncertain. Therefore, the 
difference in the faint-end slopes found between these studies shows the uncertainty involved 
in measuring background contamination when using photometric surveys in clusters (the 
optimal redshift for background subtraction is found to be z ~ 0.15 (Driver et al 1995)), 
different selection effects and methodologies. Moreover, an accurate determination of the 
Virgo LF is further complicated by the prolate shape of this cluster (Yasuda et al 1997; West 
& Blakeslee 2000; Arnaboldi et al 2002). The Ursa Major LF is found to be significantly 
different from those in richer clusters like Virgo and Coma (Trentham et al. 2001). For 
example, it does not show the steep increase at the faint-end, indicating a relatively smaller 
population of dwarfs with respect to giants. 

Recently, there have been two comprehensive studies of the field LF from SDSS (Blanton 
et al 2001) and 2dFGRS (Cross et al 2001; Norberg et al 2002). Because of the need 
for spectroscopy in estimating the LF for field galaxies, these normally span a brighter 
magnitude range than their cluster counterparts, derived from photometric surveys with 
statistical background correction. This, combined with the fact that these two studies are 
based on different filters and magnitude scales compared to present work, makes a detailed 
and quantitative comparison difficult. Nevertheless, the faint-end slope of the Coma LF 
from this study (-1.18±^) agrees with those from the SDSS (-1.24 ± 0.05) and 2dFGRS 
(—1.21 ± 0.03) field LFs. Moreover, the SDSS reveals a turn-over in the faint-end of the LF 
for red and a relatively steep LF for blue galaxies. This effect is not observed in Figure 9. 
However, in both the SDSS and 2dFGRS, strong correlations are found between luminosity 
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and surface brightness, in qualitative agreement with results from the present study. 

6. Discussion 

The shape of the LF at faint magnitudes depends on the efficiency with which gas is 
converted into stars in low mass systems (Dekel & Silk 1986; Efstathiou 2000). For example, 
dark halos in dense environments could collect gas (through intracluster pressure) and be 
turned into stars but dark halos in unevolved (diffuse) environments could not, leading to 
more low luminosity galaxies (dwarfs) per unit total mass in evolved compared to unevolved 
or field environments (Tully et al 2001; Somerville 2001). However, this effect is contrasted 
by the on-going process of galaxy interaction and tidal stripping in dense regions of clusters, 
resulting in a decrease in the number of faint galaxies and hence, flatter faint-end slope for 
core LFs, while the reverse is the case in areas with small galaxy densities (poor clusters and 
groups)- (Lopez-Cruz et al. 1997; Phillipps et al 1998; Adami et al. 2000; Conselice 2001). 
The observed faint-end slope of the LF is likely to be fixed by a competition between these 
two effects. 

The R-band LFs measured here show a slight increase in their faint-end slope from 
— 1.16^02 in Comal (inner region) to — 1.21+o;q4 in Coma2 (intermediate region) and — 1.29^ 
in Coma3 (outer region). However, this is only a la effect. Contribution from the NGC4839 
group to the faint-end slope of the Coma3 field is examined by removing galaxies associated 
with this group and estimating the Coma3 LF again. This does not significantly change the 
faint-end slope of the Coma3 field. Recently, using a spectroscopic sample of 60 clusters with 
different richness classes, selected from the 2dFGRS, De Propris et al (2003) found a com- 
posite LF with a = — 1.25±0.03, similar to that of the field LF. However, for a given spectral 
type of galaxies, they found the cluster LF to be significantly steeper than the LF for field 
galaxies of the same spectral type. The R-band luminosity of galaxies, used in the present 
study, is sensitive to the old stellar population (ie. mass function) and hence, not signifi- 
cantly affected by on-going star formation. Therefore, the similarity of the LFs between the 
core and outskirts of the cluster, and their closeness to the field LFs, estimated from SDSS 
(Blanthon et al 2001) and 2dFGRS (Cross et al 2001; Norberg et al 2002), suggests that the 
average luminosity produced for a given galaxy mass is little affected by the processes which 
make significant changes to the morphologies and spectral types of individual galaxies. This 
is in qualitative agreement with results in De Propris et al (2003). 

The similarity of the faint-end slopes of the R-band LF between core and outskirts of the 
Coma cluster is in agreement with the results in Trentham (1998) who uses a significantly 
deeper photometric survey, extending to M R = —12. However, Beijersbergen et al (2002) 
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found an increase in the faint-end slope of U-band LF with radius from the center of the Coma 
cluster, indicating an abundance of star-forming galaxies in low density regions. Considering 
the observed LFs from these studies, the theoretical prediction of a turn over in the LF at faint 
magnitudes, due to the suppression of star formation in low mass galaxies (Efstathiou 1992; 
Chiba & Nath 1994; Thoul & Weinberg 1995) is not supported. However, due to the limited 
depth of the spectroscopic survey in this paper (Mr < —16) and uncertainties in background 
correction and normalisation in photometric surveys (Trentham 1998; Beijersbergen et al 
2002), one cannot yet constrain these models. A single component Schechter function form 
gives acceptable fits to the observed LFs over a range of 7 magnitudes (—23 < M R — 
5 log h 65 < —16) at both the core and the outskirts of the Coma cluster. However, in a 
recent study, Yagi et al (2002) showed that the composite LF from 10 nearby clusters is well 
described by a Schechter form while the LFs for individual clusters are not. Furthermore, 
they found the composite LF to have a steeper faint-end slope compared to the general field, 
in contrast to the present study and the results from De Propris et al (2003). However, this 
is < 2a effect and not statistically significant. 

The characteristic magnitudes of LFs derived here agree with recent estimates from 
Beijersbergen et al (M* = -21.91+°;*j + 5 log h 65 ; M* = -20.03+°:4o + 5 lo 9 after 
converting to the waveband used in this study. However, comparison with a large sample 
of nearby (z < 0.11)- (De Propris et al 2003) and intermediate redshift (z ~ 0.3)- (Yagi et 
al 2002; Valloto et al 1997; Gaidos 1997) clusters show that while values are consistent 
(M* = -22.04 ± 0.11 + 5 log h 65 (Gaidos (1997) and M* = -22.24 ± 0.20 + 5 log h 65 
(Yagi et al (2002)), our M* B estimate for Coma is fainter by almost one magnitude (M B = 
-20.80 ± 0.10 + 5 log h 65 (De Propris et al (2003) and Valloto et al (1997))). Differences 
between the faint-end slope of the LFs between these studies have only a small effect in 
this comparison. This indicates that galaxies in the Coma are, on average, redder than 
other clusters. This is unlikely to be due to a redshift effect, as De Propris et al sample 
covers a range z < 0.11. Also, the richness of the Coma cluster cannot be responsible 
for this result, since there are also many rich clusters in the Yagi et al sample. Moreover, 
fitting a subset of their clusters with a < 800 Km/sec, De Propris et al still estimate 
M B = —20.72 ±0.10 + 5 log /i 65 , similar to that from their full sample. This is likely due to 
the very large fraction of early-type galaxies in the Coma cluster, as also shown in Figure 9. 

A dip is detected in the total LFs for the Coma cluster in both R and B-bands. This 
appears at the same magnitude, scaled by the mean B — R color of galaxies in Coma. A 
similar effect is also found in the LFs from other studies of the Coma (Biviano et al 1995; 
Driver et al 1994), other nearby (Wilson et al 1997) and intermediate redshift clusters (Yagi 
et al 2002; Dahlen et al 2003) but not in deep field surveys (Ellis et al 1996; Cowie et al 
1997). We also find the amplitude of the dip to increase from R to B-band. The presence 
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of the dip implies a two-component shape for the cluster LF, consisting of a Gausssian and 
Schechter function at the bright and faint magnitudes respectively (Biviano et al 1995; Yagi 
et al 2002). Considering the total LF in B — R color intervals in Figure 9, we find a distinct 
change in the shape of the LF for red galaxies (B — R > 1.35) at Mr = — 18 mag., both 
at the core (Comal) and outskirt (Coma2 and Coma3) fields. This result is in qualitative 
agreement with the finding by Yagi et al (2002) who attribute the change in the LF shape 
to contribution from early-type galaxies, defined as those with r 1//4 -law surface brightness 
profile. 

The question of particular interest here concerns the nature of galaxies contributing to 
the faint-end of the Coma LF and if these are different between the Comal and Coma3 
fields. Study of the LF in color intervals (Fig. 9) shows an abundance of low luminosity 
red galaxies (B — R > 1.3) in Coma 1, compared to star-forming objects (B — R < 1.3), 
with the difference decreasing towards the outskirts (Coma2 and Coma3). Moreover, the 
faint-end slope of the red component of the LF is significantly steeper (a = —1.6) than 
that observed for blue galaxies (a = —1.1). This implies that low luminosity red galaxies 
dominate the faint-end slope of the LF, in agreement with the result in Conselice (2001) 
who finds a steep faint-end slope for the LF in A0426 due to the presence of a low-mass 
red population. Therefore, the main contribution to the faint-end of the Coma LF (and in 
clusters studied by De Propris et al (2003)), is the large number of faint red galaxies that are 
present (presumably produced from the equally large number of faint blue galaxies in the 
field). A relatively higher surface density for this population in the denser region of the cluster 
(Comal) compared to the less dense region (Coma3) indicates that dynamical stripping of 
high mass systems in cluster environments is, at least partly, responsible for the formation 
of low-luminosity red galaxies. However, other mechanisms must be at work in less dense 
regions of Coma to produce this population. Evidence for gas deficiency in dwarf galaxies 
in Coma3 comes from the absence of HI emission around N4839 group (Bravo-Alfaro 2000). 
This indicates that low luminosity galaxies in the outskirts of the Coma cluster have gone 
through a gas-loss process, either through SNe winds (due to their small potential well) or by 
dynamical effects after passing through the cluster. Indeed, in a study of radial dependence 
of spectroscopic line indices of these galaxies (after correcting for luminosity dependence), 
Carter et al (2002) find a significant gredient in Mg 2 , in the sense that galaxies in the core 
have stronger Mg2 indices and hence, higher metallicities. This is caused by trapping of 
the material in galaxies in dense regions due to external pressure by intracluster medium, 
leading to higher metallicities. Moreover, X-ray morphology of the N4839 group in Coma3 
field indicates that this group is in the process of infall, showing evidence of interaction with 
the cluster (Neumann et al 2001). Further observational evidence for dynamical stripping 
in dense regions of clusters comes from the discovery of diffuse arcs at the core of the Coma 
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cluster (Trentham & Mobasher 1998; West & Gregg 1999). 

There appears to be a clear separation in luminosity distribution of galaxies as a function 
of their effective surface brightness. This monotonic surface brightness-luminosity relation 
is useful in studying properties of dwarf galaxies with respect to giants (Ferguson & Binggeli 
1994). Spectroscopically confirmed members of the Coma cluster have a relatively brighter 
effective surface brightness (/i e // < 20 mag./arcsec 2 ) compared to field galaxies, implying 
destruction of low surface brightness galaxies in environments with higher local densities. 
This is predicted by the harassment scenario of galaxy formation, that due to small potential 
well of low luminosity/ low surface brightness galaxies, up to 90% of their stellar content 
could be ejected in dense regions of clusters due to interaction with larger galaxies (Moore 
et al 1996). However, this same process also predicts formation of low surface brightness 
galaxies through stripping at the cores of rich clusters (Moore et al 1999). Therefore, the 
present data cannot be used to constrain models for the formation of low surface brightness 
galaxies. High resolution spectroscopy, measuring features diagnostic of age and metallicity 
is found to be more effective in studying formation and evolution of low surface brightness 
galaxies (Carter et al 2002). 



7. Summary 

The spectroscopic survey of galaxies in the Coma cluster, performed by Mobasher et al 
(2002), is compiled with a brighter sample (Edwards et al 2002), providing a wide-area (1 
deg 2 ) and deep (R ~ 19.5) survey, covering both the core and outskirts of the cluster. The 
final survey consists of a total of 1191 galaxies of which, 760 galaxies are spectroscopically 
confirmed members of the Coma cluster. After correcting the spectroscopic sample for 
incompletness, the LFs are constructed, spanning the range —23 < Mr— 5 log < —16, and 
covering both the core and outskirts of the cluster. Results from this study are summarised 
below: 

• The R-band LFs are similar at the core and outskirts of the Coma cluster, with no 
evidence of a steep faint-end slope, found in previous studies (mostly based on pho- 
tometric surveys). However, due to spectroscopic limitations, the current sample only 
extends to Mr = — 16 mag. while the observed up-turn in the LF slope is seen at some- 
what fainter magnitudes in photometric surveys. The total R-band LF for the Coma 
cluster, fitted to a Schechter form, is; M* R = -21.79l£{jg + 5 log h 65 ; a = -1.18lg;^. 
This is found to be close to that in the general field. 

• As the R-band mostly measures contributions from the old stellar population, rela- 
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tively unaffected by star formation, the similarity of the LFs implies that the average 
luminosity for a given galaxy mass is little affected by the processes which dictate their 
morphologies and spectral types. 

• The total B-band LF, fitted to a Schechter form, is: M* B = -19.951^ + 5 log h 65 ; 
a = — 0.96^q;q2- This shows a dip at M B = — 18 mag. A similar feature, although 
with smaller amplitude, is observed in the total R-band LF at the same luminosity 
(shifted by mean B — R color). This is likely to be due to contribution from luminous 
early-type galaxies (M B < —18 mag.) to the total LF. This suggests that the total 
LF could best be fitted by 2-components, a gaussian (at bright magnitudes) and a 
power-law (at fainter magnitudes). 

• A study of the LFs in B — R color intervals show a steep faint-end slope for red 
(B — R > 1.35) galaxies, both at the core and the outskirts of the cluster. This 
population of low luminosity red galaxies has a higher surface density than the blue 
(B — R < 1.35) star-forming population (both at the core and outskirts of the cluster), 
and dominates the faint-end of Coma cluster LF. 

• A monotonic correlation is found between the effective surface brightness and luminos- 
ity of galaxies. Cluster galaxies are found to have a higher surface brightness than their 
field counterparts. This implies destruction of low surface brightness galaxies in dense 
regions of clusters, indicating effect of local environment on formation of galaxies. 
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